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The unique properties of graphene have raised high expectations regarding its application in carbon-based 
nanoscale devices that could complement or replace traditional silicon technology. This gave rise to the vast 
amount of researches on how to fabricate high-quality graphene and graphene nanocomposites that is 
currently going on. Here we show that graphene can be successfully integrated with the established 
metal-silicide technology. Starting from thin monocrystalline films of nickel, cobalt and iron, we were able 
to form metal silicides of high quality with a variety of stoichiometries under a Chemical Vapor Deposition 
grown graphene layer. These graphene-capped silicides are reliably protected against oxidation and can 
cover a wide range of electronic materials/device applications. Most importantly, the coupling between the 
graphene layer and the silicides is rather weak and the properties of quasi-freestanding graphene are widely 
preserved. 

For almost ten years, graphene, a two-dimensional crystal of carbon atoms packed in a honeycomb structure, 
has been in the focus of intensive research. Its remarkable and promising electronic properties have raised 
high expectations regarding its potential use in the next generation electronics of a post-silicon era 1 . 
However, a graphene-based technology first and foremost requires reliable and low- cost approaches to fabrica- 
tion of high-quality and large-scale graphene layers on a variety of functional substrates. A major challenge is the 
integration of graphene into the current Si-based technology, for which original solutions have become available 
recently 2 " 6 . One basic idea is to use silicon to create a semiconducting 7 " 9 or an insulating 10 buffer between 
graphene and a metallic substrate, and thus to obtain structures that are suitable for electronic applications. 

We propose to use monocrystalline films of nickel, cobalt and iron as the substrates, on top of which high- 
quality graphene can be easily created using the Chemical Vapor Deposition (CVD) approach 11 " 15 . The silicides 
of these three metals in particular exhibit a variety of functional properties 1617 and are already extensively used. 
It is worth noting that metal silicides are generally currently employed for many applications including 
complementary metal oxide semiconductor (MOS) devices 1819 , thin film coatings 20 , photovoltaics 21 ' 22 , thermo- 
electrics 23 , catalysts 24 and biomolecule sensing 25 . 

For instance, the use of Ni silicides as metal gate electrodes in complementary MOS devices has attracted 
particular attention 26 . Ni-rich silicides, like Ni 3 Si, Ni 2 Si and Ni 31 Si 12 , can be used for p-MOS gates due to their 
higher work function, while NiSi, NiSi 2 and Ni 3 Si 2 are attractive for n-MOS applications 27 . Ni and Co silicides also 
exhibit intriguing electronic properties, especially in the form of nano wires or nanobelts 28 " 31 . The compound 
CoSi 2 , similarly to Ni silicides, is a material of choice to contact the source, drain and gate areas in complementary 
MOS 1819 . The quite different electronic and magnetic properties of Fe-based silicides 32 ' 33 ensure huge potential of 
these materials for various applications, too 34 ' 35 . For example, oc-FeSi 2 exhibits metallic properties while P-FeSi 2 is 
a narrow-gap semiconductor. Ferromagnetic Fe 3 Si layers can be used as a source for injection of spin-polarized 
electrons into a semiconductor 36 ' 37 . 

The aim of this work is to uncover the basic properties of the graphene/silicide interfaces. Here, we demonstrate 
that intercalation of silicon atoms underneath the graphene top layer and a series of annealing cycles allow one to 
create a variety of metal silicides with different properties, all nicely covered and protected by graphene. The 
graphene layer itself barely interacts with the silicides underneath and shows quasi-freestanding character. We 
anticipate a wide range of opportunities where the properties of such silicides, the graphene layer on top of the 
semiconducting silicides or of both combined could be utilized. 
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Figure 1 | Schematic illustration of the proposed approach. At first, graphene is produced by CVD on 10 nm-thick metallic films of Ni( 111), Co(0001) 
or Fe(110) grown epitaxially on W(110). Then, cycles of silicon deposition and annealing are performed to promote intercalation of silicon atoms 
underneath graphene and initiate silicide formation. Desirable composition and properties of graphene-protected silicides depend on the amount of 
deposited silicon as well as temperature and duration of annealing. 



Results 

Fig. 1 schematically illustrates the proposed approach. In the begin- 
ning, a high-quality and large-area graphene layer is CVD-grown on 
a single-crystalline thin Ni, Co or Fe film substrate. Then, silicon 
atoms are gradually deposited on top of the freshly prepared gra- 
phene layer. Their penetration underneath graphene occurs already 
at room temperature (RT). The mechanism of such intercalation 
process is based on the well-known fact that graphene grown on 
transition metal surfaces like Ni, Co, Pt, Ru etc. possesses a domain 
structure. Thus, the intercalation of different elements like Si or Au 
occurs via the domain boundaries 38 . Since the solubility of Si in Ni is 
relatively high and the diffusion barrier for the deposited atoms on 
graphene is evidently quite low, silicon atoms will migrate to these 



domain boundaries and penetrate underneath graphene already at 
room temperature. In order to promote the chemical reaction 
between Si and metal, and obtain the desirable silicide phase, sub- 
sequent annealing cycles must be applied. The stoichiometry of the 
metal silicides can be controlled by the quantity of intercalated sil- 
icon as well as the duration and temperature of the annealing pro- 
cedure. The silicon intercalation also considerably attenuates the 
strong hybridization between the electronic states of graphene and 
the substrate and, in particular cases, leads to the stabilization of the 
properties of quasi-freestanding graphene. 

Graphene on Ni-silicide structures. We first go through the detailed 
analyses of the steps towards silicide formation in the case of Ni. 



c 

CD 



c 
o 
'to 

(f) 

"e 

CD 
O 

O 




T 



287 286 



i i i i I i i i i I i i i i 
285 284 283 



i i I i i i i I i i 
101 100 



I I I I I 
99 



Binding energy (eV) 



98 



Figure 2 | Revealing the early stages of Si intercalation and silicide initiation on Ni substrate. The sequence of C Is and Si 2p core level spectra for 
graphene/Ni(lll) (a), and with increasing doses of Si deposited on top, accompanied by annealing treatments (b-f). The highlighted part of the C Is 
spectra is associated with decoupled graphene regions. 1 A corresponds to the Si coverage of 2.88- 10 14 at/ cm 2 . The C Is and Si 2p spectra were taken with 
470 eV and 135 eV photons, respectively. 



SCIENTIFIC REPORTS | 3 : 2168 | DOI: 1 0.1 038/srep021 68 



2 



Fig. 2 shows the core-level photo electron (PE) spectra of graphene 
freshly prepared on Ni(l 1 1), after deposition of small doses of silicon 
atoms, and after the post- annealing treatment. The C Is spectrum of 
native graphene/Ni at the bottom of the left panel looks very similar 
to that previously reported for this system 13 . Even a small dose (~6 
A) of silicon deposited at RT notably changes the spectral shape and 
gives rise to a shoulder at lower binding energies (BE). The Si 2p 
spectrum does not reveal a well-resolved spin-orbit doublet. Instead, 
we see a rather extended and unshaped bump feature at the BE of 
—99.3 eV, which is consistent with the reported values for Si 2p 3/2 of 
clean silicon (which is 98.9-99.6 eV depending on the surface 
reconstruction type). Together with a blurred LEED pattern (not 
shown) in comparison to that obtained from native graphene on 
nickel, these observations suggest that the major amount of the 
deposited silicon atoms creates disordered, amorphous islands on 
top. However, the shoulder in the C Is spectrum at 284.5 eV 
implies that the penetration of silicon underneath graphene 
already starts at RT, and graphene is already partly lifted up and 
detached from the substrate 39 . 

Annealing of this structure at 450°C for 40 minutes almost restores 
the shape of the C Is spectrum. At the same time, it leads to complete 
disappearance of the bump feature and the formation of a well- 
resolved Si 2p spin-orbit doublet. The BE of the Si 2p 3/2 component 
is 99.5 eV, which is in a good agreement with the known values for 
nickel silicides 40 . Compared to the spectrum before annealing the Si 
2p intensity became notably smaller and the LEED pattern again 
reveals a sharp (1X1) hexagonal structure with low background. 



Hence, the post-annealed state may be described as complete inter- 
calation of silicon atoms underneath graphene. One may assume that 
due to their notable diffusion and relatively high solubility in nickel, 
the majority of Si atoms penetrate deep into the bulk of the substrate, 
thus becoming undetectable by XPS. This could reasonably explain 
the drop in the PE intensity of the Si 2p spectrum. The restored shape 
of the C Is spectrum suggests that the chemical interaction of gra- 
phene with the underlying silicon-containing nickel substrate is sim- 
ilar to what was found for the native graphene/nickel system. 
Essentially, the intercalation of silicon atoms under graphene and 
their diffusion deeper into the Ni film leads to the formation of a Ni- 
Si solid solution in the subsurface region with a Si concentration up 
to the solubility limit of about 10 at.% 41 . Further increase of the Si 
concentration should give rise to Ni silicide formation starting from 
the Ni 3 Si phase. Thus, the Ni silicide formation should be steered by 
the balance of the intercalation rate and the bulk diffusivity. 

Upon depositing further portions of silicon at RT and subsequent 
post- annealing processes, the system strictly follows the scenario 
described above. The Si 2p spectrum taken immediately after the Si 
deposition reveals again a bump-like structure showing that plenty of 
amorphous silicon islands are formed on top of the graphene layer. 
Its PE intensity decreases considerably after the annealing cycles, and 
gradually, the well-resolved Si 2p doublet becomes dominant. The 
latter, apparently, originates from the ordered phases of intercalated 
Si atoms in nickel. This time longer annealing at higher temperature, 
500°C, does not completely remove the bump in the Si 2p spectrum. 
This may be explained by the formation of oxidized silicon islands 
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Figure 3 | Monitoring the variety of graphene-capped Ni-silicide phases. Left panel: the sequence of C Is and Si 2p spectra for freshly prepared graphene/ 
Ni(lll) (a), and after gradual intercalation of Si promoted by annealing at 430°C for 10 min for each step (b-e). The highlighted part of the C Is 
spectra corresponds to decoupled graphene regions. Right panel: deconvolution of the Si 2p spectrum (c) recorded with 165 eV and 420 eV photons. The 
surface (s) and bulk (b) components are highlighted. Si/Ni ratio is estimated from Si 2p and Ni 3p spectra and used for evaluation of the silicide 
stoichiometry. 
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Figure 4 | Evidence of electronic graphene-substrate decoupling upon 
silicides formation. NEXAFS spectra at the carbon K-edge for graphene/ 
Ni( 111) (a) , and for the systems with intercalated Si (b-d) . The part of the 
spectrum that does not exhibit angular dependence is highlighted in red. 
The geometry of the experiment is depicted in the inset. 

which were not intercalated under the given conditions. Clearly, the 
substrate becomes richer in silicon, implying the formation of the 
silicide phases underneath graphene. The LEED picture taken after 
the last treatment, corresponding to the upper PE spectra in Fig. 2, 
now reveals a sharp (2 X 2) structure with respect to the Ni(lll) 
surface. The stoichiometry of the formed silicide will be discussed 
below. 

The evolution of the system with increasing deposited, and sub- 
sequently intercalated silicon, up to 34 A of nominal thickness, is 
shown in Fig. 3. The series was measured on a freshly prepared 
graphene/nickel substrate. Each annealing step was done for 
10 min and at a temperature of 430°C. 

This time, the shoulder at the lower BE in the C Is spectrum 
shaped into a well-defined peak at 284.5 eV. It is worth noting that 
we do not detect any sign of a PE signal at —283 eV BE which would 
be indicative of carbide phases 13 , i.e. degradation of the graphene 
layer. It is interesting to compare the C Is spectra depicted in 
Fig. 2e and Fig. 3b, as the amount of deposited silicon in these two 
experiments was the same but the annealing temperatures are differ- 
ent. Apparently, reducing the temperature as well as the duration of 
annealing effectively increased the amount of silicon in the subsur- 
face region underneath graphene. This is reflected in the C Is spec- 
trum as a noticeable enhancement of the PE intensity at 284.5 eV. 
The LEED picture again reveals a sharp (2 X 2) pattern. Taking into 
account the silicon concentration estimated from the Si 2^/Ni 3p 
intensity ratio we suppose, that such structure corresponds to the 
formation of ordered Ni 3 Si silicide under graphene. The solid state 
diffusion and the formation of a bulk alloy rather than a surface 
compound is supported by the fact that Si concentration decreases 
after annealing. This process is quite similar to the recently observed 
formation of the bulk Ni 3 Al phase with a (2 X 2) surface structure as 
a result of graphene intercalation with Al atoms 42 ' 43 . The possible 
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Figure 5 | Quasi-freestanding character of graphene on silicides. The 

sequence of ARPES data for graphene on Ni(l 1 1) (a), and at its detachment 
upon Si intercalation and silicide formation (b-d). The data were obtained 
near the K-point of graphene Brillouin zone using 40 eV photons. 

structural model is shown in Fig. SI. The Si 2p spectral shape 
becomes very different from that of the Si-depleted systems described 
above. We see now a new well-defined peak at lower BE. To gain 
insight into the spectral structure of the discussed Si 2p core-level 
spectrum, we recorded it at two different photon energies, 165 eV 
and 420 eV, and then deconvolved the obtained Si 2p spectra as 
shown in the right panel of the Fig. 3. 

Both PE spectra were normalized to the maximum intensity. Each 
spectrum can be fitted with two doublets, where the two components 
likely reflect electron emission from Si atoms in the subsurface and in 
the bulk. The lower BE component is enhanced at the lower photon 
energy, which allows it to be assigned to the surface structure. Hence, 
the higher BE component originates from silicon atoms embedded in 
the bulk of the system. 

We analyzed the Si 2^/Ni 3p intensity ratio at each stage in order to 
evaluate the stoichiometry of the formed structure and to follow its 
evolution during the experiment. The obtained data together with 
the well-known phase diagrams for the Ni-Si system suggest that 
starting from the solid solution we have passed from the Ni-rich 
phase Ni 3 Si to NiSi. We have to admit that there may be a certain 
gradient of the silicon concentration from the surface to the bulk. 
Thus, for each particular system we may have several Ni- silicide 
phases 44 with higher Si concentration at the near-surface region, 
which is probed with XPS and solely defines the properties of the 
graphene-silicide interface. It should be noted that at low dosage of 
the intercalated silicon the multicomponent structure of the C Is 
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Figure 6 | Formation of graphene-capped silicides in the Co-based system monitored with XPS. Evolution of the C Is and Si 2p spectra with Si 
intercalation underneath graphene on Co(0001) (left panel). The highlighted part of the C Is spectra corresponds to decoupled graphene regions. 
Separation of the silicon atoms embedded in the cobalt layer near the surface and in the bulk obtained via deconvolution of the Si 2p spectrum (right panel). 



spectrum indicates coexistence of the graphene/silicide and gra- 
phene/nickel interfaces. 

In Fig. 4, we show the evolution of NEXAFS C-K edge spectra, 
which allows us to follow possible modifications of the n and a 
electronic states of graphene with increasing intercalated Si. The 
shape of the C-K edge spectrum of freshly prepared graphene on 
Ni(lll) has been intensively discussed recently 45 . It contains two 
characteristic features, a double peak at —285 eV and a manifold 
structure at —292 eV. They reflect C Is core electron excitations into 
the 7i* and a* states, respectively. Both features are substantially 
broadened compared to what is observed for quasi- freestanding gra- 
phene because of the strong hybridization between the electron states 
of carbon and nickel. The red- colored parts of the spectra do not 
show the strong angular dependence typical for the n* and a* fea- 
tures. This suggests that the sp 2 character of the n states is altered due 
to the interaction with the Ni atoms below. Upon the transition from 
the solid solution via Ni-rich phases to the NiSi silicide, the lineshape 
of the spectrum becomes more and more similar to that of the free- 
standing graphene. This suggests that with increasing amount of 
silicon atoms inside nickel, the Ni 3d states participate more actively 
in the formation of the Ni-Si bonds in silicides 46 and, in turn, interact 
less and less with graphene. As a result, the electronic coupling 
between the graphene and the silicide substrate becomes weak. 

This is clearly seen in Fig. 5, where we show the changes in 
the electronic band structure of graphene/Ni(lll) with increasing 
intercalated Si. The band structure for the freshly prepared gra- 
phene/Ni(lll) is already well-known and has been a subject of 
extended discussions 47 " 49 . Due to the strong chemical interaction with 
the substrate, the n band of graphene is strongly shifted in energy. 
The Dirac-point lies about 2.6 eV below the Fermi level. However, 



even a small amount of intercalated Si, corresponding to the Ni-rich 
silicides, strongly changes the topology of the electron bands. The n 
band is shifted back to lower BE again, and the cone apex appears 
quite close to the Fermi level. At the same time, the PE intensity of the 
respective Ni 3d states considerably attenuates. With further increase 
of the intercalated Si and stabilization of silicide phases, the Dirac- 
point reaches the Fermi level, i.e. its position in the freestanding 
graphene. 

Graphene on Co- and Fe-silicide structures. Having established a 
clear picture of the formation of Ni-silicides under CVD-grown 
graphene, let us now turn to a brief description regarding Co- and 
Fe- based systems. 

We performed the experiments in the same way and deposited 
first the silicon atoms on top of the freshly prepared graphene layer 
on monocrystalline films of cobalt or iron at RT. Later on, the anneal- 
ing cycles were applied and the full process was monitored by ana- 
lyzing the C Is and Si 2p spectra, which are depicted in Fig. 6 and 
Fig. 7, respectively. The stoichiometry of the metal- Si systems was 
evaluated via the Si 2p/Co(Fe) 3p ratio. 

Discussion 

Looking closely at the presented spectra, we may see and discuss 
certain similarities as well as differences of the silicide formation in 
these three systems. We should note that the silicide formation is 
solely defined by the local subsurface Si concentration, which, in 
turn, is governed by the balance of intercalation and diffusion rates. 
The local concentration of silicon atoms should be higher than that 
defined by the silicide stoichiometry, for example, not less than 50 
at.% for NiSi. Both diffusion and intercalation rates increase with 
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Figure 7 | Formation of graphene- capped silicides in the Fe-based system monitored with XPS. Evolution of the C Is and Si 2p spectra with Si 
intercalation underneath graphene on a Fe(llO) layer (left panel). The highlighted part of the C Is spectra corresponds to decoupled graphene regions. 
Separation of the silicon atoms embedded in the iron layer near the surface and in the bulk obtained via deconvolution of the Si 2p spectrum (right panel). 



temperature, these dependencies are described by Arrhenius equa- 
tions with two different activation energies. The choice of the heat 
treatment temperature and duration can be performed using these 
equations. Unfortunately, the lack of fundamental information on 
this matter allows us to explain the behavior of these three systems on 
the phenomenological level only. 

Despite the fact that Ni and Co metals exhibit very similar prop- 
erties, the difference in the behavior of graphene/Ni and graphene/Co 
systems upon Si intercalation can be explained as follows: the Ni 3 Si 
phase exists in a wide temperature range in the Ni-rich part of the Ni- 
Si system. In contrast, the Co 3 Si phase is observed only within a short 
range of relatively high temperatures. The Ni 3 Si phase has a similar 
crystal structure to Ni, i.e. it arises due to the ordering in Ni-Si solid 
solution. In such systems typically the interactions between different 
atoms are stronger than for the atoms of the same element. This 
enhanced interaction, in turn, promotes interdiffusion due to ther- 
modynamic reasons 50 . Because of this, the Si diffusion in Ni is notably 
faster than that in Co in the corresponding metal-rich solid solutions. 
Therefore, in the Co-based systems the surface- related component of 
the Si 2p spectrum always has higher intensity than the bulk one, 
which is not the case for the Ni-based structures. Note that for Fe- 
systems the situation becomes somewhat intermediate, because sim- 
ilar to Ni, Fe forms extended solid solutions with Si up to the ordered 
intermetallic phases Fe 3 Si. Hence, the discussion given above is true 
for the graphene on iron system, too. However, in contrast to nickel, 
iron has a higher melting point. Therefore, the bulk diffusion is acti- 
vated at essentially higher temperature. Thus, at the same temper- 
ature, Si evidently has higher diffusivity in nickel crystal than in iron 
one. This is in line with our spectral observations described above. 



In summary, we have demonstrated that the ordered silicides of 
Ni, Co and Fe with different and well-defined stoichiometries can be 
produced under CVD -grown graphene. The obtained materials are 
reliably protected by the graphene layer, making them compatible for 
a variety of nanoscaled applications at ambient conditions. Our find- 
ing suggests that the graphene layer is weakly electronically coupled 
to the silicide phases, thus keeping its remarkable properties ready 
for use. We anticipate that the presented results will create a solid 
platform to further explore the opportunity to integrate graphene 
with the mature silicon-based technology. We believe that this would 
give rise for further exploitation of the unique properties of different 
metal silicides under graphene's umbrella as well as the properties of 
graphene itself for the widely discussed applications in the next- 
generation electronic devices. 

Methods 

The studied graphene samples were synthesized in situ by CVD from propylene 
(C 3 H 6 ) on the Ni(lll), Co(0001) and Fe(llO) surfaces. For this purpose we prepared 
metallic films with a thickness of ~ 10 nm grown epitaxially on the W( 110) substrate 
in ultra-high vacuum conditions. The high crystalline quality of the studied samples 
was confirmed by the sharp low energy electron diffraction (LEED) patterns. After 
that, graphene layer was grown on nickel and cobalt surfaces at a pressure of pro- 
pylene of ~ 1 • 10" 6 mbar and the substrate temperature of 600°C during 15 min. At 
these conditions the reaction is well-known to be self-limited to single graphene layer. 
The parameters of CVD for the iron-based system were slightly increased to avoid 
carbidization and produce a complete graphene layer, namely 5-10" 6 mbar, 700°C 
and 25 min, respectively (see Supplementary Note 1 for details). The penetration of 
the silicon atoms into the metal substrates and the subsequent formation of various 
types of silicides underneath the graphene top layer were studied using X-ray 
photoelectron spectroscopy (XPS), angle resolved photoelectron spectroscopy 
(ARPES), near-edge X-ray absorption fine structure spectroscopy (NEXAFS), and 
LEED. The well-known stability of graphene-based systems to ambient atmospheric 
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conditions 51 " 53 was proved for the studied structures, too. Experiments were per- 
formed at BESSY synchrotron radiation facility - at the beamlines U125/2-SGM, 
UE56/2-PGM-2 and at the Russian- German beamline. ARPES measurements were 
performed at room temperature with a Scienta R4000 spectrometer. 
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